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Effect of a variable cosmological constant 
on black hole quasinormal modes 
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Many different cosmological models have been proposed to address the cosmological constant 
problem and the coincidence problem. We compare here four different models that can be used 
to describe an effective (time-dependent) cosmological constant A(z). A numerical analysis of the 
A(z) evolution obtained for each model shows that it can be used for distinguishing between all four 
models. We calculate next the <u(A) frequencies for quasinormal modes of gravitational perturbations 
of Schwarzschild-de Sitter black holes at different redshifts. Considering that the variation of A 
happens on cosmological timescales, the combined <u(A(z)) could be used in principle to track the 
evolution of the cosmological constant. We quantify the resulting minute frequency shift in the 
quasinormal mode frequencies and show that the relative frequency shift grows as M 2 . However, 
even in a most optimistic scenario with an extremely high mass supermassive black hole there is no 
prospect for the detection of this effect. 

PACS numbers: 04.70.Bw,95.36.+x,04.25.Nx 


I. INTRODUCTION 

According to recent observational evidence, the uni¬ 
verse is spatially flat and is currently undergoing an ac¬ 
celerated expansion phase 01- The simplest physical 
explanation for this result would be the existence of a 
cosmological constant driving the current acceleration of 
the universe, but there is a huge difference between the 
observed value of the cosmological constant and the vac¬ 
uum energy density derived by quantum field theories 
(the “cosmological constant problem”). Another prob¬ 
lem with the cosmological constant model is the “coin¬ 
cidence problem”: why should the accelerated expansion 
occur exactly now in the history of the universe? 

These two problems have motivated many different 
cosmological models, ranging from scalar field models of 
dark energy to coupled dark energy models and mod¬ 
ified gravity alternatives [§]. So many different mod¬ 
els, based on different physical assumptions, must be 
constrained by observations. Usual tests include obser¬ 
vational constraints from type la supernova data, cos¬ 
mic microwave background radiation and baryon oscilla¬ 
tions. However, the comparison between different models 
should be done carefully, otherwise the results could be 
model- or calibration-dependent Q. 

An effective cosmological “constant” can be derived 
from dynamical dark energy models, its value varying 
slowly over cosmological timescales. For physical pro¬ 
cesses happening at shorter timescales, this variation will 
be negligible. We study here how the quasinormal modes 
of gravitational perturbations of a Schwarzschild-de Sit¬ 
ter (SdS) black hole will be affected by a variable effective 
cosmological constant A as a function of the redshift z. 
The quasinormal mode frequencies will be a function only 
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on the parameters of the spacetime 0 [§]. In the timescale 
for quasinormal mode oscillations, A will be essentially a 
constant, but black holes at different redshifts will have 
different quasinormal mode frequencies. 

In Section mi we present four different cosmological 
models and numerical results for each A (z). The vari¬ 
ation of the SdS quasinormal modes as a function of A is 
shown in Section IIIII In Section HVl we propose a cosmo¬ 
logical application of our results, calculate the resulting 
frequency shift for supermassive black holes and discuss 
the detectability of this effect. Finally, we give our final 
remarks in Section IVl 

II. COSMOLOGICAL MODELS 

We consider here four different cosmological models 
(models A-D), that have different underlying physical as¬ 
sumptions. From each one we can derive equations that 
describe the cosmological constant as a function of the 
redshift. We present numerical results for each model in 
Fig- HI For all models we keep z < 15, which we consider 
as an upper bound for black hole formation [7|. 

We do not intend for this Section to be a comprehen¬ 
sive review on these models, so we will point the reader 
to the original papers in each subsection for more details. 

A. Dimensional analysis 

In |8j], a model is proposed following some simple as¬ 
sumptions on quantum cosmology and a dimensional 
analysis. In this scenario, the effective cosmological con¬ 
stant is a function of the scale factor a as A oc a~ 2 or, as 
an explicit function of the cosmological redshift z, 

A(z) = A 0 (l + z) 2 , (1) 

where Ao is the present value of the cosmological con¬ 
stant. This would lead to an early universe value of A 
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several orders of magnitude larger than the present Ao, 
which is a useful feature in light of the cosmological con¬ 
stant problem. This model is presented as “model A” in 
Fig. HI 


B. Coupled quintessence 


The coupled quintessence model proposed in jo] has a 
light scalar field <f> (the quintessence) with an exponen¬ 
tial potential U(<p) = Ae^ 2 ^^^, linearly coupled to the 
ordinary matter, that is responsible for the cosmic ac¬ 
celeration as a dynamical cosmological constant. This 
model is written in terms of the variables 


K (f> 



na I p r 



( 2 ) 


where H is the Hubble parameter and p r is the radia¬ 
tion density. By using the variables xi, X 2 and £ 3 , the 
Klein-Gordon equation for the scalar field and the Fried¬ 
man equation for the spacetime can be re-written as the 
following system 

x'l = X 1 ^ - lj - px\ + /3(1 -x\-xl~ £ 3 ), 

( x ' 

Xn = IIX 1 X 2 +£ 2 ( 2-1 -- 

V *3 

^3 =-y(l-3x?+ 3 xl-x\), (3) 

where the independent variable is a = log a, the prime 
denotes d/da and /3 is a dimensionless constant given by 
/? = 1 / 3 / 2 C/K (C is the coupling constant between the 
scalar field and the ordinary matter). 

We remark here that the parameters /3 and p, (related 
to the scalar field potential) are all that is needed to 
completely specify the model. However, the numerical 
solution found is strongly dependent on the initial values 
supplied for ([3]). The initial values (at early times) re¬ 
quired for cosmic solutions are near (xi, £ 2 , £ 3 ) = ( 0 , 0 , 1 ) 
and the exact values used here were found by trial and 
error, until we approached the standard ACDM values 
for the present cosmological constant density parameter, 
Ha 0 ~ 0.7, and for the present matter density parame¬ 
ter, f2 mo ss 0.3. With these values, the model makes the 
transition from a radiation dominated universe, at early 
times, to a matter dominated universe until close to the 
present time, when we have a A-dominated stage. 

With the numerical solution of the system Q, we can 
then obtain the scalar field density parameter as a 
function of a , 



= x\ + x\ = — , (4) 

Pc 

where p\ = A/8irG and the critical density p c can be 
written as 

Pc = ^( l + z)\ (5) 

x 3 


where p ro is the present radiation density. Therefore, we 
can explicitly write 


A(z) = A 0 (l + z) 4 


x i 


x\ 


l 3 0 


( '2 0 


( 6 ) 


with A 0 = 8nGp ro (x? + x 2 0 ) / x 2 0 and the Xj 0 are the 
present values of the Xi at 2 = 0. This model is presented 
as “model B” in Fig. [L] 

Since there is no analytical solution for A (z) in eq. 
we performed a 4th order polynomial fit for the numerical 
results, 


—r-^ = b 0 + 6 i(l + z) + 62(1 + z) 2 + 

Ao 

+ ^ 3(1 + ~) 3 + 64(1 + z) 4 , (7) 

with the coefficients given in the left side of Tabic [I] 


C. Interacting holographic dark energy model 


In 0 , a holographic dark energy model with an inter¬ 
action between dark energy and dark matter is studied. 
The future event horizon is chosen as the infrared cutoff, 
providing a time-dependent ratio of the matter energy 
density and the dark energy density. This model also al¬ 
lows the transition of the dark energvequation of state 
to phantom regimes, as proposed by [llj . 

For this model we can write from the Friedman equa¬ 
tion, 


Wd 

% 

ir_ 


(1 — f Id) 


fip 
1 2-^/TTd 

3 + 3c 


+ 



3 b 2 
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Qd 
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( 8 ) 

(9) 


where fin is the dark energy density parameter, c is a 
constant and b 2 is the coupling constant for the interac¬ 
tion between the dark energy and the dark matter. After 
obtaining numerical solutions for fljj and H , we can fi¬ 
nally write 


A = A 0 


h 2 n D 
HZShTo' 


( 10 ) 


where A 0 = 377qHd 0 and H 0 is the present value of 
the Hubble parameter H(z). This model is presented 
as “model C” in Fig. |T] 

Following the procedure we used in III B1 we did a poly¬ 
nomial fit of the numerical results obtained, 

= C 0 + Ci(l + z) + C 2 (l + z ) 2 + 

Ao 

+ c 3 (l + ~) 3 + £ 4(1 + z) 4 , ( 11 ) 

with the coefficients given in the right side of Table [I] 
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bo 

1.2161(58) 

co 

0.6726(18) 

bi 

- 0.2763(49) 

Cl 

0.1936(15) 


0.0800(12) 

C2 

0.11372(39) 

63 

0.03661(12) 

C3 

- 0.002830(37) 

64 

-5.964(36) x 10~ 4 

C4 

6.76(12) x 10" 5 


TABLE I. Values for the 6; and a coefficients for the model B 
and model C fits given in eq. 0 and eq. GD. respectively, and 
their standard deviations, calculated with a standard least 
squares fit using the numerical solutions presented in Fig. 0 


values grows to ~ 15% at 2 = 15. For z 3.5, model 
B is intermediate between models C,D and model A, be¬ 
coming larger than models C,D by a factor 2 at z w 7. 

Based on these results, we can conclude that all mod¬ 
els are distinguishable from each other solely by their 
A(.z) predictions. This could lead to a very direct way 
of comparing cosmological models, provided we have ob¬ 
servational data on A at different redshifts. In the next 
Section we will discuss how black hole oscillations could 
help us to obtain these results. 


D. Cosmological particle creation 

Dark energy is represented by a dynamical cosmologi¬ 
cal constant in [l2j . This model proposes that the space- 
time expansion process can extract non-relativistic dark 
matter particles from the vacuum, at a constant rate T 
given by 


T = ^(1 — ^m o )H 0 , (12) 

where Q mo is the present matter density parameter. As 
a consequence, it can be shown that A = 2 TH and we 
can write it explicitly as 

A(z) = A 0 [l — Q mo + fi mo (l + z) 3 /-]. (13) 

where A 0 = 3(1 — V, mo )HQ. The concordance value ob¬ 
tained for f i mg in this model with a joint analysis of the 
matter power spectrum, the position of the first peak of 
the CMB anisotropy spectrum and the Hubble diagram 
for type la supernovae is fl mo ~ 0.45 (higher than the 
usual ACDM value). 

This is another example of a model with interaction in 
the dark sector. We can also cite as a motivation the fact 
that particle creation is expected to happen in curved 
spacetimes. The coincidence problem is considered to 
be alleviated in this context, because of the consequent 
slower decay of the a ratio. This model is presented 

as “model D” in Fig. [T] 


E. Comparison between models A-D 

The solutions obtained for A(z) in the models A-D 
described in sections III AHII Dl are presented in Fig. |T| 
We can see that, in model A, A grows much faster than 
the others, being larger by a factor 2 already at z « 1 . 
For 0 < z 2.5, the magnitude of A follows model B < 
model C < model D. For 2.5 ;$ z 3.5, we have model 
C < model B < model D. Finally, for z 3.5 we have 
model C < model D < model B. 

Models C and D, even though motivated by different 
physical assumptions, provide very close results for the 
entire z range considered. The difference between their 
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FIG. 1. A/Ao as a function of z for models A-D discussed in 
sections iHAHirDl The analytical solutions for models A and 
D are given in 0 and G3- Model B is given by eq. m, 
where we used the numerical solution of the system m with 
/r = 0.1 and /? = 0.5. Model C is given by eq. 03, where we 
used the numerical solutions of eq. 0 and eq. © with c = 1 
and b 2 = 0.04. 


III. QUASINORMAL MODES FOR THE 
SCHWARZSCHILD-DE SITTER BLACK HOLE 


There is a vast literature on the quasinormal modes of 
stars and black holes, see |5, 6 ; , 13] for some very good 
reviews on the subject. The problem of finding the quasi¬ 
normal modes of the SdS black hole has already been 
previously considered in the literature, see for instance 
PlUS!, an d it has gained more attention recently in the 
context of black hole thermodynamics and area quanti¬ 
zation (see 0,0 and references therein). We review 
here only very briefly the most important points needed 
for our analysis. 

The metric of a Schwarzschild-de Sitter spacetime with 
black hole mass M and cosmological constant A is de¬ 
scribed by 


ds 2 = -f(r)dt 2 + fir)- 1 dr 2 + r 2 dil 2 , 
where dfl 2 = dd 2 + sin 2 9dtp 2 and /(r) is given by 


,, , , 2 M r 2 

fi r ) = 1 ~ A— 


(14) 

(15) 



















4 


The radial component of a perturbation in this back¬ 
ground satisfies a wave equation that can be put in the 
following form 

^ + (‘^ + F ' (r) )' , ’ , = 0 ' <16) 

where x is the tortoise coordinate given by dr/dx = f(r ) 
and the form of the potential Ve(r) depends on the nature 
of the linear perturbations considered (see, for instance, 
the discussion in 0). 

The quasinormal modes are solutions to the pertur¬ 
bation equation m that satisfy the boundary condi¬ 
tions corresponding to purely outgoing waves at infin¬ 
ity and purely ingoing waves at the black hole horizon. 
These solutions are characterized by complex frequencies 
to = uj r + iuii (the quasinormal frequencies) that depend 
only on the parameters of the background spacetime, and 
not on the details of the initial perturbation. 

For gravitational perturbations in the SdS spacetime, 
we can fit the dependence of the dimensionless quasinor¬ 
mal frequency Mu> as a simple polynomial function of 
A: 


Ad lo — o >q -(- mA -f- ( 12 -^ U3A 3 T a^A 1 ^ , (17) 

where the complex coefficients a, are given in Table HU 
Our empirical fits for u> r and u>i are presented in fig. [2j 
together with numerical results taken from 0. We chose 
a quartic fit based on the behavior of the data, as can be 
seen in the figure. We do not propose this empirical fit 
with a deeper analytical justification, but as a useful tool 
for our following analysis. 


a,i 

Re (a;) 

Im(ai) 

a o 

0.3730(45) 

-0.0888(12) 

ai 

-0.87(70) 

0.13(18) 

02 

-53(28) 

14.0(75) 

a 3 

9.0(40) x 10 2 

-2.3(10) x 10 2 

a4 

-5.3(18) x 10 3 

1.39(46) x 10 3 


TABLE II. Values for the a,i coefficients for the empirical fit 
given in eq. 03 and their standard deviations, calculated 
with a standard least squares fit using the data from 0 
shown in fig. [2] present in c = G = 1 units. 


Combining the empirical fits (El) with the solutions we 
obtained in Section [TT] for A(z), we can produce uj r (A(z)) 
and Wi(A(z)). These combined expressions show how 
the quasinormal mode frequencies would be affected by 
a variable cosmological constant: black holes at differ¬ 
ent redshifts will have different corrections from A(^), 
that will depend on the particular cosmological model 
used. The necessary assumption implicit for this analy¬ 
sis is that the timescale for the evolution of A is much 
longer than that of the black hole perturbations, which 
allows us to use the metric m with the variable A(^). 
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FIG. 2. Upper plot: Muj r as a function of A The data points 
are taken from [Ijj and the solid line represents the empirical 
fit II171 1. with the values for the coefficients given in Table 1X11 
Lower plot: same as the upper plot, but for Muji. 


IV. A COSMOLOGICAL APPLICATION 


Since the actual value of Ao is so small, we will keep 
now only the terms up to first order in A in the absolute 
frequency difference caused by the time dependence of 
the cosmological constant, that can be written as 

Auj r = (u r - u 0r ) = Re(ai)^ - 1^ , (18) 


where to o r = tu r (Ao). If we return now to physical units, 
we can assess exactly how small this effect is and the 
feasibility (or not) of using it as a real cosmological test. 
For A = 0, the frequency / measured in Hz is 


r _ Re(ao) / c 3 \ 
2t r \GMJ 


12.05 



kHz, 


(19) 


and restoring the physical dimensions to A and Re(ai) 
gives 
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A = A {4m) = A * 10_7m_2 ' (2 °) 

/ M \ 1 2 * * K o 

Re(oi) = Re(ai) — 2.18 x 10 6 m 2 , (21) 

\ m qJ 


(recall that a±A is dimensionless). If we take now the cur¬ 
rent value of the cosmological constant as approximately 
Ao = 10 _ 52 m~ 2 , and ai = —0.87 from Table ITT1 then we 
have aiAo = —1.9 x 10and we can hnally 
rewrite eq. m as 


A/ = -6.14 x 10 “ 42 



A(s) 

Ao 


1 Hz. 


( 22 ) 


which allows us to write the relative frequency shift 


A/ 

I 


= -5.10 x 10 “ 46 



Ao 



(23) 


Equations E^l) and (12311 confirm the intuition that the 
effect of the cosmological constant on the quasinormal 
modes should be very small, but they go beyond this 
intuition and quantify the effect 0 We can see that it 
becomes larger at larger masses, but even in a case with 
a M = 10 1o Mq black hole and an enhancement factor 
of 100 from the evolution of A there is no prospect of 
measuring A/. In this case, we would have the relative 
frequency shift A/// a: 10” 24 . 


V. FINAL REMARKS 

We have presented here a comparison between four 
different cosmological models with variable cosmologi¬ 
cal constant. We presented analytical solutions for A(z) 
when possible, in the case of models A and D (this re¬ 
sult is new and was not presented in 0) and numerical 
solutions with corresponding polynomial fits for models 
B and C, which are also novel results not provided in 

references [9j and S3- 

Our results show that all four models considered here 
can be distinguished from each other by their predicted 
evolutions, as can be seen in Fig. |T] Even for models C 
and D, which show the closest behavior, this difference 
amounts to ~ 15% at z = 15. For the other models, the 
difference can be as large as a factor ss 10 at the same 
redshift. 

We also studied the influence of the variable cosmolog¬ 
ical constant in the quasinormal modes of gravitational 
perturbations in a Schwarzscliild-de Sitter black hole, and 
presented a numerical fit for w(A). An increase in the 
value of A will decrease both uj r (and increase the oscil¬ 
lation period) and |wj| (and increase the damping time 
of the perturbation). 

The resulting effect in the quasinormal mode frequen¬ 
cies is extremely small. We quantified this effect and cal¬ 
culated it for a supermassive black hole as a cosmological 
application of our results. However, the detectability of 
the effect of the evolution of the cosmological constant in 
black hole quasinormal modes is very far from the current 
status of any gravitational wave detectors. 
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